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Abstract 

Through an Hamiltonian action we write down the system of equations of motions 
for a mixture of thermocapillary fluids under the assumption that the internal en- 
ergy is a function not only of the gradient of the densities but also of the gradient of 
the entropies of each component. A Lagrangian associated with the kinetic energy 
and the internal energy allows to obtain the equations of momentum for each com- 
ponent and for the barycentric motion of the mixture. We obtain also the balance of 
energy and we prove that the equations are compatible with the second law of ther- 
modynamics. Though the system is of parabolic type, we prove that there exist two 
tangential acceleration waves that characterize the interfacial motion. The depen- 
dence of the internal energy of the entropy gradients is mandatory for the existence 
of this kind of waves. The differential system is non-linear but the waves propagate 
without distortion due to the fact that they are linearly degenerate (exceptional 
waves) . 
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1 Introduction 



Liquid-vapor and two-phase interfaces are generally represented by a material 
surface endowed with an energy related to Laplace's surface tension. In fluid 
mechanics and thermodynamics, the interface appears as a surface separat- 
ing two media. This surface has its own characteristic behavior and energy 
properties [Ij. Theoretical and experimental detailed studies show that, when 
working far from critical conditions, the capillary layer has a thickness equiv- 
alent to a few molecular beams [2J. 

Molecular models such as those used in kinetic theory of gas lead in interfacial 
layers to laws of state associated with non-convex internal energies, e.g., the 
van der Waals models [3~Pf5] . These models appear advantageous as they pro- 
vide as even more precise verification of Maxwell's rule applied to isothermal 
phase transition [BJ. Nonetheless, they present two disadvantages: 
First, for densities that lie between phase densities, the pressure may become 
negative. Simple physical experiments can be used, however, to cause traction 
that leads to these negative pressure values [7118] . 

Second, in the field between bulks, internal energy cannot be represented by a 
convex surface associated with the variation of densities and entropy. This fact 
seems to contradict the existence of steady equilibrium state of the matter in 
this type of region. 

To overcome these disadvantages, the thermodynamic investigative replaces 
the non-convex portion corresponding to internal energy with a plane domain. 
The fluid can no longer be considered as a continuous medium. The interface 
is represented as a material surface with a null thickness. In this case, the 
only possible representation of the dynamic behavior of the interface is one of 
a discontinuous surface, and its essential structure remains unknown. 

In the equilibrium state it is possible to eliminate the above disadvantages 
by appropriately modifying the stress tensor of the capillary layer, which is 
expressed in an anisotropic form. As a consequence, the energy of the con- 
tinuous medium must change [5 p| 9]. A representation of the energy near the 
critical point therefore allows the study of interfaces of non-molecular size. 
This approach is not new and, in fact, dates back to van der Walls [10] and 
Korteweg [TT]; it corresponds to what is known as the Landau-Ginzburg the- 
ory [5]. The representation proposed in the present study is based on the 
notion of internal energy which is more convenient to use when the tempera- 
ture is not uniform. One of the problems that complicates this study of phase 
transformation dynamics is the apparent contradiction between Korteweg's 
classical stress theory and the Clausius-Duhem inequality [12J. Proposal made 
by Eglit [13], Dunn and Serrin [H], Casal and Gouin [15] and others rectify 
this apparent anomaly for liquid-vapor interfaces of a pure fluid. 

To study capillary layers and bulk phases, the simplest model in continuum 
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mechanics considers an internal energy e as the sum of two terms: a first one 
corresponding to a medium with a uniform composition equal to the local one 
and a second one associated with the non- uniformity of the fluid [3"]nU] . The 
second term is approximated by a gradient expansion, typically truncated to 
the second order. In the simplest version of the theory we have e = pa(p, s) + 
| m ( Vp) 2 , where p is the matter density, s the specific entropy, a the specific 
internal energy of the fluid assumed to be homogeneous and m is a coefficient 
independent of s, Vp and of any higher derivatives. Obviously, the model is 
simpler than models associated with the renormalization-group theory |16j . 
Nevertheless, it has the advantage of easily extending well-known results for 
equilibrium cases to the dynamics of interfaces [T7lfT8] . In such a model, s 
varied with p through the interface in the same way as in the bulk, then s 
would always be that function of p which, at given temperature T satisfied 



With this assumption, s = s(p) and e = e(p), so that the points representing 
single-phase states in the p, s, e space lie on a curve instead of on a surface 
e = e(p,s). This was the original assumption of van der Waals which was 
later justified by Ornstein in 1909 for a system composed of molecules with 
long-ranged but weak attractive force; this assumption is not exact for more 
realistic potentials^}. 

As coexistence curves these are no way peculiar; the only peculiarity is that the 
single-phase states -in this version of this approximation- have collapsed onto 
the coexistence curve instead being represented by the points of an extended 
two-dimensional region of which the coexistence curve is merely a boundary. 
There is then no proper two-density description of the one-phase states of a 
one-component system in the lowest order of the mean-field approximation. 
There is such a description of the two-states, where not even in mean-field ap- 
proximation is there any discernible peculiarity; but in practice the potential 
for the two-density form of the van der Waals theory is then not constructed 
by the prescription in (pQ) but by other means. For example Rowlinson and 
Widom introduce in [I], chapter 9, an energy arising from the mean-field the- 
ory and depending on the deviations of the densities s and p, say, from their 
values at the critical point and the gradients of these densities. It is also seen 
that in c-components systems, c + 1 densities -the densities p 1; ...,p c of the c 
components and the entropy density s- may vary independently through the 
interface. 

Aside from the question of accuracy, there are also qualitative features of some 
interfaces in physical-chemistry, especially in systems of more than one com- 
ponent, that require two or more independently varying densities for their de- 
scription. An example is strong positive or negative adsorption of a component 

1 Main sentences in this paragraph and more precisely the comments on static 
interfaces are issued from the book of Rowlinson and Widom [3] and its bibliography. 
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i associated with a non-monotonic profile p^iz) where z is the spatial variable. 
In the one-density theory based on the approximation ^ — T = and = 
for all j 7^ i, the resulting one-density model leads inevitably to a monotonic 
Pi(z). In a theory based on two or more densities, by contrast, we may have 
a realistic trajectory with which is associated non-monotonic behavior of one 
or more of the components if we suppose Pj(z) to be monotonic. 
We must also allow the independent variation of at least two densities, entropy 
included, if we are to account contact angles in three-phase equilibrium : J. 
Cahn made the remark that we might use a two- (or more-) density van der 
Waals theory to describe the case of non-spreading in the three phase equilib- 
rium, ([3], chapter 8). Then, at equilibrium, Rowlinson and Widom pointed 
out that for single fluids the model must be extended by taking into account 
not only the strong variations of matter density through the interfacial layer 
but also the strong variations of entropy. 

Also in dynamics, for an extended Cahn and Hilliard fluid, the volume inter- 
nal energy e is proposed with a gradient expansion depending not only on 
grad p but also on grad s corresponding to a strong heat supply in the change 
of phases: e = f(p, s, Vp, Vs). The medium is then called a thermocapillary 
fluid [1320]. 

To extend the model to fluid mixtures corresponding to c-component systems 
and realistic potentials in molecular theory of fluid interfaces, the internal en- 
ergy is assumed to be a functional of the different densities of the mixture. 
In all the cases where strong gradient of densities occurs - for example shocks 
or capillary layers - the internal energy is chosen as a function of successive 
derivatives of densities of matter and entropies. To be in accordance with the 
physical phenomena presented in [I], one will consider the internal energy of 
a two-component fluid mixture with an internal energy depending also on the 
gradients of entropy of each component. The internal energy is a Galilean in- 
variant, it does not depend on the reference frame; hence the internal energy 
depends also on the relative velocity between the two components of the mix- 
ture. 

The conservative motions of thermocapillary fluid mixture are relevant to the 
so-called second gradient theory [21] and we obtain a complete set of bal- 
ance equations for conservative motions; we extend this result to the dissi- 
pative case. Our goal is simply to verify the consistency of our model with 
the second law of thermodynamics. We consider a special case of dissipative 
thermocapillary mixtures where the introduction of dissipative forces is only 
done in the framework of the first gradient theory; we deduce the Fick law 
as a consequence of a friction behavior between the components and from the 
equations of motion of the components. In such a case, extended thermody- 
namic principle (as Gibbs identity) provides a set of equations that satisfy the 
entropy principle, thereby making these irreversible motions compatible with 
the second law of thermodynamics. 
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The idea of studying interface motions as localized traveling waves in a multi- 
gradient theory is not new and can be traced throughout many problems of 
condensed matter and phase-transition physics [22] • In Calm and Hilliard's 
model [3J, the direction of solitary waves was along the gradient of density 
[TTf22] . The introduction of the model of thermocapillary fluid mixture pro- 
vides a better understanding of the behavior of motions in fluid mixture inter- 
faces: it is possible to obtain the previous solitary waves but also a new kind 
of adiabatic waves may be forecasted. These waves are associated with the 
spatial second derivatives of entropy and matter densities. For this new kind 
of adiabatic waves, the direction of propagation is normal to the gradient of 
densities. In the case of a thick interface, the waves are tangential to the inter- 
face and the wave velocities depend on the constitutive equations. Finally we 
observe that, also if the differential system associated with the wave motions 
is non-linear, the waves propagate without distortion due to the fact that they 
are linearly degenerate (exceptional waves) [23] • 



2 Equations of thermocapillary mixtures 

2. 1 Conservative motions of thermocapillary mixtures 

To derive the governing equations and boundary conditions in the dissipative- 
free case, we use the Hamilton principle of least action |24| . In continuum 
mechanics, the principle with a system endowed with an infinite number 
of degrees of freedom was initiated by Lin [25], Herivel [26], Serrin |27j . 
Berdichevsky [2H] and many others; it was proposed by Gouin for fluid mix- 
tures [21] . The main idea is to propose a Lagrangian which yields the behavior 
of the medium as the difference between a kinetic and a potential energy per 
unit volume. Then, the variations of the Hamilton action obtained as a lin- 
ear functional of virtual displacements allow to find the governing equations 
and boundary conditions. For real media, the irreversibility is introduced both 
in equations of motion and equation of energy by using a classical approach 
through the dissipative function, diffusion and heat fluxes. 

We study a mixture of two fluids: the motion of a two-fluid continuum can be 
represented using two surjective mappings (i=l,2) 

(t,x) ^X. = Ai(t,x) 

where (t, x) belongs to [ti,t 2 ] X D t , a set in the time-space occupied by the 
fluid between times t\ and t 2 . Variables Xj denote the positions of each com- 
ponent of the mixture in reference spaces D 0i . Variations of particle motions 
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are deduced from virtual motions 



where scalars Xj are defined in a neighborhood of zero; they are associated 
with a two parameter family of virtual motions of the mixture. The real motion 
corresponds to x, = 0, the associated virtual displacements are 



l*s=o 



They generalize what is obtained for a single fluid [30]. To the virtual dis- 
placements <5jXj, we associated its image in the physical space D t occupied 
by the fluid mixture at time t [24T3T] . 

<9x 

-SiXi. 



dX, 

Conservation of matter for each component requires that 

p i det(J^=p io (X i ), (2) 

where p io is the reference volume mass in Dot and det f^fr) t ne Jacobian 
determinant of the motion of component i. In different iable cases eqs (|2J) are 
equivalent to the equations of balance of matter densities p i 

^ + divfou,) = 0, (3) 

where Uj denotes the velocity vectors of each component i. Now, we assume 
that the mixture has an entropy for each component [32]; for conservative 
motions, the equations of conservation of specific entropies Sj are 

^ + div(p^u,) = 0. (4) 

Then, relations 

define an isentropic motion of the fluid mixture. We deduce the following 
relations of tensorial quantities [21129] 



6iUl = ^d±~ dx ^' 6iPi = ~ diV 6iSi = ~dx^' ^ 

di d 

where — = — + u,.V denotes the material derivative relatively to the corn- 
eft dt 

ponent i. We assume that the volume potential energy of the mixture is in the 
form 

e = e{p i ,s i , Vpi, Vsi,w), 
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with w = ui — U2 the relative velocity of the two components of the mixture. 
This means that the fluid mixture is a function not only of the densities of 
matter p i and specific entropies Sj but also of the gradients of p i and Sj. The 
fact that e depends on two entropies is classically adopted in the literature 
[29|33|34j . Moreover, for a two- velocity medium, there is no coordinate sys- 
tem within the framework of which any motion could be disregarded. So, the 
standard definition of potential energy leads to its dependence on the relative 
motion of the components. The dependence of e with respect to the relative 
velocity is analog to take into account the added mass effect in heterogeneous 
two- fluid theory as done by berdichevsky [2B] and Geurst [2335] • Let us note 
we can assume also that e is depending on (t, x); by this way, we introduce 
directly the extraneous potential of the body forces. 

The potential e is related with the volume internal energy w of the mixture 
through the transformation 

de 

W = £ - — W, 

aw 

so that, 

2 1 

e = E nPi U i + W 

2 = 1 Z 

is the total energy of the system pHf37] . 

The equation of motion of component % is given by a variational method 
associated with a Hamilton action; the vector field x 6 D t — > C li is two time 
continuously differentiable. The Lagrangian of the mixture is 

2 1 

i=l z 

and consequently, the Hamilton action between the times t\ and £2 is 

-t 2 



L dxdt. 

h JD t 

From the definition of virtual motions, we obtain immediately two variations 
of the action of Hamilton associated with i = 1,2, 



^•jSi ^i^i €,Si^$iSi,'yJ dx.dt , 

where subscript 7 corresponds to spatial derivatives associated with gradient 
terms; as usually summation is made on repeated subscript 7 from 1 to 3; 
1 / <9e \ ^ 

Kj = Uj + (— 1)* — 7^ — where index T denotes the transposition. Then, by 

Pi V 5w / 
integration by part we obtain, 
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* r m , 



5iI = It Id \\2 Ui ~ e ' Pl + ( e ' p ^ rt ) 5iPi + PiKl ~< 5iUi ~< 

- (e, Si -(e, SiiT ) i7 ) diSi - {e tP . Jip % ) tl - {e, ai JiSi) tl ) dxdt. 

Let us denote by 

„ de be . N 

PA = and ^ = ( 6 ) 

dsj dpi 

where <9 is the variational derivative operator. That is to say, 
P^i = e, Sl -(e, Si , 7 ) l7 = e, Si - div and h { = e tP . - (e il0 . J >7 = e, p . - div 
with, 

\I>, = — and <fr, 



Introducing Ri = \u\ — hi, and taking into account of the expressions for 9i 



dVsi dV Pl 

— h- • 

and hi given by eqs (jBJ), we get 

8J = ( [ {Ri 5iPi + p^Kj • 5jUi - pfii 5iSi - div (<&j 5^ + *i ^s^)} dxcft, 
and from relations (jSJ), we obtain 



+ P#i -g^ d - div (*i 5iPi + *i 5iSi) > dxdt. 
Consequently, 



t2 f { ( dRi : n dsi di~Kj T ^ T dui 



(7 1 

+ ^(^Kf Ci) - div (p^ Ci - Pi^i Kf Ci + hPi + *i ] d*dt. 
The Stokes formula and relation ([3]) yield 



Jti JdDt 

( do ■ ds ' i 

-n. f Pi R % Ci - PiUi Kj Ci-^i^Ci- *i ^ Ci - Pi^ div C 4 J ^x^, 
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where dD t (of mesure da x ) is the boundary of D t , n is the unit external normal 
vector to dD t and g is the velocity of dD t . If we consider a vector field x G 
D t — > and its first derivatives vanishing simultaneously on the boundary 
o\D t , the Hamilton principle expressed in the form: V 5ja = leads to 

and consequently, 

Let us note that the value of the first member of eq. (|S)) is equal to 
and eq. (jHJ) yields 



= p^j V Si - pj V /ij . 

Taking into account of eq. ([3]) of conservation of mass of component i, we 
deduce 



^ + dlv(ftU , Su , ) + ( _ ir ( dlv ,(*)%|(*) T 



and due to the fact that 

divui f-^-J u i = div f f ^- J ®Ui,. 

y ow y ax y aw y I y aw y 

we get the equations of motion of the two components in the form 
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h (Ih [ ( J^-j ®U, I I = p,6>, V.s, -ftV/!, (9) 



We consider only the isotropic case where the potential energy e of the mixture 
can be written in terms of the isotropic invariants 

Pij = VprVp^ Xij = Vp r Vs i? 7 - = Vsi-Vsj, (i, j = 1,2) and u = ^w 2 . 

£ = ^{.Pii s i) Piji Xij i liji k"0> 

Then, the equation of motion of each component of the mixture is 



(), '< Ui + <liv(p,u, u,) • (10) 



dt 



d , x ( dui\ 



T 



- -1 )' ( ^(ow) + a ^— J w + div(a w (g) U;) | /.</•', V«, - /.., V //,. 
where a = — . 

In this case, \I>i and <frj can be written 

2 2 

= E A, v Pi + / -o v - s / , = E v Pj + vs j , 

3=1 3=1 

with Cij = (1 + Sij) e,p.. , Dij = e, x .. , E^ = (1 + <%) e, 7tf , where is the 
Kronecker symbol. 

The simplest model is when Cy = Cjj, Dy, E 1 ^- = i?^ are constant. Then, 

2 1 1 
e = e( Pi , Si , w) + g ^ V ^ ' V ^ + D « V ^< ' Vs i + 2 ^' Vs * ' VSj (U) 

where the associated quadratic form with respect to the vectors Vpj and 
Vsj is in the form Y^j=i \ C i j^p i ■ Vp^ + D^ Vpj • Vsj + | E^ Vsj • Vsj. 
This quadratic form is assumed positive such as the effect of gradient terms 
increases the value of the internal energy with respect to a mixture in a ho- 
mogeneous configuration. 
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2.2 Equation of total momentum and equation of energy for conservative 
motions of thermocapillary mixtures 



We limit first to the conservative case. We notice that the equation of motion 
of each component is not in divergence form. Nevertheless, by summing eqs 
(II 01) with respect to i, we obtain the balance equation for the total momentum 
in a divergence form. In fact, eqs ffTUl) imply 



fe + div (Pi u i ® u i)) - div ( a w ® w) = 

a (£) W+ (EpAVs.-^V^ (12) 
In coordinates, the second member of eq. ( fT2l) is 

2 

a w u w 1)U + Y, z H s ix - (e s . J ^ s f)7 - ft (e, p .) + ft (e ft> J , 

i=l 

where v is summed from 1 to 3. Noting that 

2 

e ,7 = E e >«i S *>7 + W S VT + 6 ,ftPi,7 + e ,Pi,„Piry + a W v W l,» , ( 13 ) 

8=1 

we obtain 



■4i 

a w u w^ u + Y e - 4 Si, 7 - ( e , Si ,J „ s i,7 - P< ( e ,P.) + A ( e ,ft,„ J = e ,7 - 

i=l 

2 

E £ ,V S W + 6 ,ftPi,7 + € ,fii,*Pi,m + ( e ,*iJ „ S ^7 + A ( e ,P») , - Pi ( e ,ft,„ 



i=l 

= e n + E (-Pi e A + A ( e ,/v) J - ( $ ^Pi, 7 + ^ s i,7y .. • 

i=l V ' 7 '7 

and consequently the equation of motion for the total momentum is 
dpu 



dt 



div (p Uj (g) Uj) — pa w £§) w — cr^ = 



(14) 



where p = p x + p 2 is the total volume mass, pu = p x u.i + p 2 u 2 is the total 
momentum and cr = + cr 2 is the total stress tensor such that 

<7iry = (~Pi + Pi div $i) <^ 7 - $^p i)7 - ^s i)7 , with P, = p^ - 
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Let us notice that if e depends also on (t, x) corresponding to an external force 
potential, an additive term appears as body force in relation ([13]) and in eq. 
(1141) the body force appears in the second member. This is not the case in 
eqs ([H]) which include the body forces coming from e (depending on (t,x)) in 
terms hi. For the sake of simplicity, we do not introduce the body force in 
equation of the total momentum and equation of total energy 

The equation of energy of the total mixture is obtained in the divergence form. 
Let us define 



Mj = + div(p iUl ® Ui) + 



i=i 



e — a w 2 



div ^ fe Pi ( K i ■ U i - - U 



£ u-U 



where U = 5Zj=i ( ~df~'^ i J corres P on d s ^° the interstitial working: 

in the same way as for Cahn and Hilliard fluids, an additional term that has 
the physical dimension of a heat flux must be added to the equation of energy 



Theorem: For all motions of a thermocapillary fluid mixture, the relation 

2 ^ 



E-S + J2[~u 2 i-hi + 6 l s i ) Gi-Mi-m^O (15) 
is identicaly satisfied. 

The proof comes from the following algebraic calculation: 

Mf u, - Qu 2 - h^j d + pA^ = J^ u * ) + M(\prf)»i) + h i°i 
dsi dhi \i I f 9a w\ T T dui / T \ . 

+Pi e i-djr + P^ Ui + [~dT j + a w dx~ + v a u * w ) I u * 
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Let us note that 



| + div (j>, A ,„,) . a +-L'Jg + <"v (tafe) 

dpi ( de \ dsi de _dsi , , 9^ 

" dlv Uv^J + e "< aF + ^ v + * dlv (w) + 



yields 



^ (dpi de \ dsi de dsi dh { 



i=i 

Consequently, 



2 , ,j 



i=l 



£ ( Mf iK - ( -u? - fa + 9 iSi ) d ) + S = 



(/ 9a w\ T T dui / n 

[ — — I Uj + ow — Uj + div (a Ui\v 

^ ((^ \ PiU *) + 6 ~ a w2 ) + div ( K ^ Ui ~~ ^ u2 + ^ ) Ui 



Taking into account of the relations 



^l$. = ^$._$.^ u . and ^*. = ^v|>._v|>.^i u . 

or at dx or at ax 



and the definition of the total stress tensor, 



CT = £ -J2(Pi ^ - ft div$ 4 ) Id - - ^ 



i=l 



dx. dx 



where Id is the identity tensor, we deduce immediately the algebraic identity 

am. □ 
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COROLLARY: All conservative motions of a thermocapillary mixture satisfy 
the equation of energy balance 

s((g^ u +£ - aw2 ) + 

div ( ft (k, ■ u f - iuj - <t ( ) Ui j + e u - uj = 0. (16) 

This result comes from the simultaneity of relations Gi = 0, = and 
Mi = 0. 

Let us note that e — a w 2 corresponds to the volume internal energy w of the 
mixture. 

2. 3 Dissipative motions of thermocapillary mixtures 

The conservative fluid mixture model presented in section 2.1 is relevant to 
the so-called second gradient theory [21], [38]. In our form of equations of mass 
conservation for each component l^j, equation of the total momentum of the 
mixture (TBI) , equation of energy ([TBI) , the diffusion term J = p x (ui — u) does 
not directly appear but is deduced respectively from the velocities and densi- 
ties of the components. 

Our aim is to verify the consistency of our model with the second law of 
thermodynamics. The introduction of dissipative forces is simply done in the 
framework of the first gradient theory [2TJ : the dissipative forces applied to the 
continuous medium are divided into volume forces if and surface forces asso- 
ciated with the Cauchy stress tensor erf. Then, the virtual work of dissipative 
forces STi applied to the component i is in the form 

where <5Tj is a differential form. For such dissipative motions, no production 
of masses due to chemical reactions appears. 

For the same virtual displacement of two components, £ = ^, (i = 1,2) , 
the total virtual work of dissipative forces is 

When £ is a translation, the work 5T is equal to zero and consequently, 

E ff = or f 2 d = -if = i d (17) 

i=i 
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We specify later the behavior of forces ff (they will be associated with the 
diffusion term) and stress tensors erf. Taking into account of the dissipative 
effects, the equations of motion for each component become 



+ div(p iUi <g> xii) + (-1)' |^(aw) + a (^^J w + div(a w <g> m) 

= p i e i Vs i - Pi Vh i + divof + tf (18) 

Taking into account of relation flTjl . the equation of the total momentum 
writes 



dpu 

~dt 



+ div ^2 (p Ui (g> Uj) - pa w ® w - cr - cr^ =0 (19) 



with cr d = erf + cr d . 

The introduction of the heat flux vector q and the heat supply r comes from 
classical methods in thermodynamics [32,34,39j. If we write, 



M? = M, - div <r? - 

2 

S d = S - r + div q + u < - tr ((T d Aj), 



a=l 
2 



E d = E - r + div q - ^ div of Uj, 

a=l 

with Aj = - ( — — + f — ^ ] | represents the velocity deformation tensor of 
2 \ ax \ ax J I 

each component, relation ffl5l) writes as 

2 

£ d - S d - Y,Mf ■ Uj - (K 4 • m - i2j - <S< ) Gi = (20) 
i=i 

which can be considered as the dynamic form of the Gibbs identity. 

For the components of the mixture, equations of momenta and equations of 

masses are in the form 

M d = 0, Gi — 0. (21) 

The Gibbs identity ([20]) and eqs ([21]) imply S d = E d . If we assume that S d = 0, 
i.e., 

E (^f + div ^ u *)) e * + ^ u * ~ tr A *) - r + div q = °> ( 22 ) 

it is equivalent to write E d = 0, i.e., 



15 



div (feft ( Ki ■ Ui - \ u 



e - a w 2 + (23) 



2 _ 



<Xj — <r^ ] Uj I + e u — U + q — r = 0. 



Eq. (1221) is the equation of entropy and Eq. ( 1231) is the equation of energy. 

In the conservative case, the system is closed with two different temperatures 
6i (i = 1,2). In the dissipative case we need additional arguments to obtain 
equations for each entropy Sj which could remplace equations @. A possibility 
is to consider the case when the exchanges of momentum and energy between 
the two components are rapid enough to have a common temperature (this is 
not the case of heterogeneous mixtures where each phase may have different 
pressures and temperatures [10] )• Then, in dissipative case, if we know all the 
dissipative functions, the governing system is closed. Note also that it could be 
possible to consider a common temperature and entropy both for conservative 
and dissipative case. This case is connected with a conservative equation for 
the common entropy s (see Appendix). Another possibility is to assume that 
the entropy is transported along the ith component (say, for example i — 1) 
which was used for quantum fluids by Landau [JT] 

dpi 8 i- i \ 
—z— + div ( Pl s ui) = 0. 

In this case, the independent functions are Pl , s, u 1; p 2 , u 2 , where p i are sub- 
mitted to the constraints © and the case of Helium superfluid is a special 
case of our study corresponding to s 2 = and s — si. Using of this hypothesis 
is nevertheless doubtful for classical fluids. 

Hence, we may suppose a common temperature only for dissipative case. This 
means that 

9 = a = 6 2 (24) 
This hypothesis closes the system 



Now we focus on the governing equations for each components of the mixture: 

Mf = (2 = 1,2). 
For slow motions, we rewrite these equations in the following form : 

m*- Pi vhi- Pi ev Si - div af - if = o. 

If we consider the case when the motion of each component is regular enough, 
an approximative case is the case of solid displacements for the motion of each 
component; then, div erf ~ and consequently, 
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where \i i = hi — Osi is the chemical potential of the component i of the mixture 
at the temperature 6. 

Considering the difference M2 — Mj and using relation (TTTj) . we obtain, 

2_l E i^ = Vft 25 

P2 Pi PlP2 

where /z = /i 2 — fi 1 is the chemical potential of the mixture [12E3] . Let us 
introduce the diffusion flux J, 

J = pJm - u) = ^w. (26) 

p 

Equation ( J25l) implies, 

f d . w = Vp-J. (27) 



The term V /i ■ J corresponds to the entropy production due to the diffu- 
sion process. Expression (1271) is the connection between the mechanical drag 
force between components of the mixture and the thermodynamical process 
of diffusion; consequently, Eq. ( f2~2l yields 

£ (^T + div ^ u 0) 0< - *r (o? AO - r + div q - V • J = °- 

We have obtained the equations of balance of masses, equation of total momen- 
tum, Eq. ( f25i) between the components and equation of total energy by using 
an energetic method. They are the extension of classical mixture equations to 
equations of mixtures of fluids involving density gradients. 



2.4 Second law of thermodynamics 



Relation (122]) may be rewritten in the form 

(E^ + divCp^+divl-^i^tr (*f A0)+f.w-w£ 

This last equation yields the entropy production due to the diffusion, viscosity 
and heat flux processes. Due to 9 > 0, if we assume 

f> (<^A^+f d -w-W-|>0, (28) 

we get the Clausius-Duhem inequality in the form 
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Let us notice that (J2HD is satisfied if tr {erf A 4 ) > 0, f d -w > and V0- q < 0. 
The following constitutive laws are classically adopted for erf, f d and q : 

The stress tensor erf is a symmetric isotropic tensor function of Aj such that 
Yli=itr (erf Aj) ^ (a special case of such a function is erf = Aj(ir Aj)Id + 
2/^ A; with fa > and 3Aj + 2^ ^ 0). 

The heat flux vector satisfies the Fourier law 

q = — x V# with x ^ 0- 

The linear approximation known in the literature as the Stokes drag formula 
is adopted [4T1T42] 

f d = kw, k ^ 0, with f 2 d = -if = f d and w = ui - u 2 . 

Let us note that relations (|25|) and (|26|) together with the drag formula yield 
the property of the diffusion flux 

J = ~V//, (29) 

which is the general form of the Fick law [39ll^T] . So, the Fick law is not directly 
a linear phenomenological law as the Fourier law but a direct consequence of 
equations of motion and the Stokes drag force hypothesis which was previously 
noticed by Bowen for a different model |42j . 



3 Acceleration Waves 

Now, we will consider conservative motions of thermocapillary mixtures only. 
As it is well known, wave phenomena - in particular discontinuity waves (waves 
across the front of which some derivatives of the field variables have jumps) 
- are typical of models that are described through hyperbolic differential sys- 
tems. A classical example of discontinuity waves in continuum mechanics are 
the so-called acceleration waves in which, among the other variables, the ac- 
celeration jumps across the front while the velocity is continuous |44j . 

Dissipative systems and in particular models for diffusive processes have usual 
differential system with a parabolic structure and discontinuity waves are not 
admissible. A typical example of non admissibility is the one of Navier-Stokes- 
Fourier fluids. In this possible approach to obtain hyperbolic system is 

the method of the Extended Thermodynamic theory [32], valid also for rarefied 
gases. 



18 



Nevertheless in parabolic systems some discontinuity waves may propagate for 
particular initial data. The aim of this paper is to prove that for the present 
model of thermocapillarity fluid binary mixtures there exists the possibility of 
propagation of two tangential acceleration waves, provided that the internal 
energy is at least a function of the entropy gradient of one component. For this 
aim, we first briefly recall some very well known questions about discontinuity 
waves. 

A wave is a discontinuity wave if the wave front with Cartesian equation 
<f)(t, x) = separates the space in two subspaces in which there exists regular 
solutions of the differential system but across the normal direction of the front 
some derivatives of the field suffers a jump [23] . 
As usual, we indicate the jump with a square bracket, 



lcj>=0- 



^=o-t 



and we introduce the map between (t, x) and (0, £) , where £ = £(£, x) repre- 
sents the tangential manifold of the wave surface in time-space. Therefore the 
assumptions for the discontinuity waves are expressed for a generic function 
/ in the form, 

Qkf 



If] = 0; 



7fe 



Vfc; 



where there exists p ^ 1 such that 
for 1 < j < p — 1 ; 



(30) 



5 k f = 



Qkf 



7^ for k > p. 



Taking into account fl30|) and the Hadamard lemma 



we have 



QPf 



dx. 



7i 



• dx-. 



5 p f n h 



QPf^ 

dtp 



-xys p f 



where A and here n = (rii) are respectively the normal velocity and the unit 
normal vector to the wave front. 



The advantages of the previous symbols are that there exists a chain rule 
between the field derivative in the differential systems and the corresponding 
jump relation, 

d t — > — A 5; <9 7 — > n 7 5 . (31) 

We apply now this procedure to our differential systems assuming that across 
the wave front p { , Sj and its first derivatives are continuous and there are 
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jumps for the second derivative (p = 2) while the velocity is continuous and 
suffer a jump in the first derivative (p = 1 ) 

[Pi] = [*] = [Sp t ] = [Ssi] = 0; [8 2 Pl ] ± 0; [8 2 s % ] + 0; 

(32) 

[Ui]=0; W^O; i = l,2. 
5.i Jump conditions: 

From the balance of mass of the two components ([3]), taking into account of 
eqs (I3"2"j) and the chain rule (13"T1) . we obtain that the normal components of the 
first derivative of the velocities are continuous, 

[5u in ) = 0, (33) 

where Ui n = n • Uj. 

If we differentiate the entropy balance law of each component (jl]), with respect 
x, we obtain for the discontinuities 

- Vi5 2 Si + Vsi ■ 5ui = 0, (34) 

where 

Vi — A - U in 

are the relative velocities of the wave front with respect to the fluid compo- 
nents. From the balance of momentum we obtain 

-PiVidui + {-l) l a {vi5(u 2 - Ui) + (8\ii ■ (ui - u 2 )) n} 

-p i {Vs i [e i ]-[Vh i ]} = 0. (35) 

From eq. ([6]) we obtain 

2 

B t = - Pi [9i] = ( D H &Pj + E oi 52 s j) 
2 

A, ee -\ytn] =nJ2 ( c a 53 Pj + D a 53 s i) • ( 36 ) 
3=1 

Then relation (1351) becomes 

-p^Sui + {-l) l a {Vi5(u 2 - Ui) + (5uj • (u x - u 2 )) n} + SiVs; - p^A; = 0. 
If we multiply by n and we take into account of relation (I3"3"|) . we obtain 

PiAi ■ n = {-l) l a 5u { ■ (ui - u 2 ) + ^Vs, • n, (37) 
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then, we get the final jump conditions from the momentum equations 

PiViSui + (-l) l a^5(ui - u 2 ) - BiV t Si = 0, (38) 

where 

V tSi = V Si - {V Sl ■ n)n (39) 

denotes the tangential component of the gradient of entropy of each compo- 
nent. 

Therefore we obtain the algebraic system of 8 equations (1331) . (|34|) and (1381) 
for the 10 scalar unknowns 5 2 p i , 5 2 Si and <5uj (i = 1,2). Consequently, we 
needs two more conditions that are obtained by compatibility conditions com- 
ing from boundary conditions associated with the equation of motion for each 
component. In fact, we notice in Appendix that if p i5 Sj, Vpj, Vsj are con- 
tinuous through a surface of weak discontinuities, then div <frj must be also 
continuous through the surface 

[div $i] = 0. 

We notice additively that these two conditions are compatible with the Rankine- 
Hugoniot conditions associated to the total momentum (fl4"|) and the total 
energy balance law (|T6l) : In fact, eq. (fT4l) yields 

[<Ti + <T 2 ] = 0, 

while from eq. f|T6|) . we get 

[(TiU ln + Cr 2 U2n} = 0. 

Then, we obtain the two supplementary equations 

2 

[div *J = Yl {Cij 5 2 Pj + D i:i 5 2 s 3 ) = (40) 
i=i 

and the system for the discontinuities becomes an homogeneous closed system 
of 10 equations for 10 scalar unknowns 5 2 p i , 5 2 Si and <5iij [i = 1, 2) in the 
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form, 



[5u in ] = 

PiViSui + (-l) l a^5(ui - u 2 ) - Ej=i (Dji S 2 pj + Eji 5 2 s^ V t Sj = 

(41) 

— v i5 2 Si + Vsj • Sui = 

E? =1 (Cy 5 2 Pj + Dij 5 2 s 3 ) = 

We observe that the conditions fl36l) . (1371) are constraints for the jump of the 
third derivatives of densisties and entropies S 3 pj, S 3 Sj. 

Now we consider the weak discontinuities near the equilibrium of the fluid 
mixture; then 

U; = 0. 

Consequently v = A is the velocity of the acceleration wave. Due to the fact 
that thermocapillary mixtures can be considered as a mathematical model for 
interfacial layers between two mixture bulks [29J, the gradients of tensorial 
quantities Pi and Sj are orthogonal to the interfacial layers and consequently 
are collinear, 

Vsi = b Vs 2 . 

Then the second and the third equations of system fl4"Tj) allow to eliminate Uj 
and to get 

2 

c 2 ((a - Pl ) 5 2 Sl - a b 5 2 s 2 ) + £ 5 2 Pj + E jX 5 2 s 3 ) (V^) 2 = 

3=1 

and 

2 

c 2 s (-a b 5 2 Sl + b 2 (a - p 2 ) 5 2 s 2 ) + £ (D i2 5 2 p,- + £ j2 5%) (V <Sl ) 2 = 

j'=i 

where 



Therefore we obtain a system of compatibility between the variables 8 2 p 3 - and 
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5 2 Sj in the form 



C n 5 2 Pl + C 12 5 2 p 2 + D u 5 2 Sl + D 12 5 2 s 2 = 

C 12 5 2 p 1 + C 22 5 2 p 2 + D 21 5 2 s x + D 22 5 2 s 2 = 

< Du 5 2 Pl + D 21 5 2 p 2 + (E n + (a - Pl ) c 2 ) 5 2 Sl 

+ (E 12 -abc 2 s ) 5 2 s 2 = 

Di 2 5 2 p x + D 22 5 2 p 2 + (^12 - a 6 c 2 ) 5 2 s x 

+ (E 22 + (a-p 2 )6 2 c 2 ) 5 2 s 2 = 



(42) 



Let us denote 



C = C 



E = E 



T 



C*ll C*12 
C*12 C*22 

Du £"i 2 

Di 2 D 22 



D 



D xl D 12 
D 2 i D 22 



B = B J 



Pl — a a& 
a6 (p 2 — a)6 2 



6 2 P 



5 2 Pl 


, 5 2 s = 


5 2 Sl 


S 2 p 2 




5 2 s 2 



and system (j42p is writting 



C <5 2 p + D S 2 s = 



D T 5 2 p + (E - c 2 B) 5 2 s = 



A - c 2 B) S 2 s = 0. 



(43) 



which implies 
where 

A = E — D T C X D . 
From eq. (1431) it is simple to verify the property, 

Theorem: // EC — D 2 is positive definite and if we consider small diffusion, 
i.e. a < p* with 

P* 



P1P2 



Pl + P2 
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then, A and B are both symmetric and definite positive, all the eigenvalues c 2 s 
of eq. (l4*3|) are positive and two discontinuity waves exist. 

3.2 Exceptional waves: 

As in the hyperbolic case, a wave is exceptional or linearly degenerate (see e.g. 
W) if 

5\ = 0. (44) 

In this case the wave behavior is similar to the behavior in linear case and we 
do not get any distortion of the wave or shock formation. 

It is simple matter to prove that both the waves fulfill the exceptionality 
condition. In fact, taking into account of relations (13T I) . ( |39l) . we have 

5V t s = 

and from (1431) we obtain that A is function of the modulo of the tangential 
gradient of entropy 

A = A(|V fS |) 

and then (jS]) holds. 



4 Results and discussion: 

In this paper we prove that the model of thermocapillary fluid mixtures with 
dissipation yields a system of equations of motions compatible with the second 
law of thermodynamics at least in simple dissipative cases. The equation of 
motion and the equation of energy of the barycentric motion of the mixture 
are in a divergence form in conservative cases. 

Consequently, Hamilton's principle applied to fluid dynamics is a direct and 
systematic method to obtain the equations of conservative motions. This prin- 
ciple extended to each component of a mixture of conservative fluids is able to 
deduce the same number of balance equations than unknown functions. The 
method yields a non ambiguous framework for the case of non- conservative 
mixtures (with viscosity, diffusion and heat transfer). Non additional assump- 
tion but constitutive behavior compatible with the second law of thermody- 
namics is necessary. One obtains the dynamic Gibbs relation and Fick's law 
as a consequence of governing equations. 

We have seen that the dependance of an entropy gradient is necessary for the 
existence of isentropic waves of acceleration along the interfaces: the fact that 
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the internal energy depends not only on the gradient of matter densities but 
also on the gradient of entropy, yields a new kind of waves which does not ap- 
pear in simpler models. They are exceptional waves in the sense of Boillat and 
Lax [23] and they appear only in, at least, systems with two dimensions. These 
second order waves are of weak energy and consequently they are not easy to 
show up. Recent experiments in space laboratories in micro-gravity conditions, 
for carbonic dioxide near its critical point, have showed the possibility of such 
waves [16]. The experimental evidence of such adiabatic waves with other 
physical reasons we have presented in the introduction should strengthen the 
necessity to take into account of the dependence of entropy gradients together 
with density gradients in the expression of the internal energy for continuum 
models of capillarity and phase transitions. 
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5 Appendix: 



5. 1 Compatibility conditions for weak discontinuities 



We obtain also the compatibility due to the boundary conditions. For this aim 
we rewrite the variations of the Hamilton action when the equations of the 
motion (jSJ) are verified. Then, from relation ([7|), we obtain, 



= r i g Pt Kj c, 

Jh JdD t 



dpi 



ds,. 



-n. Pi Ri d - Kf Ci ~ *< -w 1 d ~ *» ~a Cj Pi®i drv d da x dt 



dx. 



dx 



25 



For a vector field x £ D t — > Ci vanishing with its first derivatives on the 
boundary dD t , we deduce immediately on a surface of discontinuity S t (where 
p i , Si and its first derivatives are continuous and there are jumps for the second 
derivative) the value of the variation of the Hamilton action, 



5J 



t-2 



9 



n. 



PiRi d - Pi^i Kf C. - ®i -q£ d - *i d - Pi®i div Ci 



da T dt 



= - f I n . \piRi] d da x dt. 
Jti JSt 



Due to the fact that 5jJ = for a vector field x £ D t — > Ci? we obtain 
[Pi-Ri] = 0, and deduce the compatibility conditions (jjOj) across the wave front, 



[div * f ] = £ (Cii S 2 Pj + A, 5 2 s 3 ) = 0. 

3=1 



5.2 Mixture with an entropy for the total fluid 



The equations of balance of matter densities p, p 1 are in the form 

dp 

-£ + div(pu) = 0, 



^ + div(p lUl ) = 0, (45) 

where p and u denote the density and the velocity of the total mixture, p x 
and ui denote the density and the velocity of one of the two components (it 
is equivalent to consider the total density p and the concentration c = p x j p 
betwen the two components of the mixture such that eq. (145!) is equivalent 
to d(pc)/dt + div(pcu) = 0). For conservative motions, the equation of 
conservation of the total specific entropy s is 

dps 

— h div(psu) = 0. 

We consider a volume potential energy of the mixture in the form 



£ = e(p, Pi, s, si, Vp, Vp 1: , Vs, u). 
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with uj = — w 2 , where w = u — u x . As in the section 2, it is possible to deduce 
the equations of motions through the Hamilton action. We denote 

. de , de , de de 
pO = , h = — , hi = — — and a = ——. 
ds dp dp 1 du 

(In fact, the potential energy of the mixture depends only on the insentropic 
invariants (Vp) 2 , (Vp^ 2 , Vp • Vp 1; (Vs) 2 ). By analogous calculations asso- 
ciated with a Lagrangian of the mixture in the form 



T 1 2 

and virtual motions such as X = 3(t, x,x) and Xi = S 1 (t,x,x 1 ), we obtain 
as in section 2, the equations of motion 



dpu . . , , 

+ div pu ® u) + 46 

at 

— |-f (aw) + a^— w + div(a w ® u) ) = p9 V s — pV /i 
ax J 

and 

— (aw) + a— - w + div(a w ® m) + p x V /i x = (47) 
at ax y 

with 

p6> = pe iS — div /i = pe iP — div $ , and hi = Pi£, Pl — div $i , 
where 

cte de de 

* = o , * = o , *1 



<9Vs aVp ' dVpi 

By summing eqs (I46f47p . we obtain the balance equation for the total momen- 
tum in a divergence form 

dpu 

— — + div (p u ® u — pa w ® w — er) = 

where cr is the stress tensor such that 



0V 7 = (-P + p div ®+pi div $x) c^ 7 - $„p - <&i v pi - ^ v s 



7' 



where P = e — pe iP — Pi£ ^ • The equation of energy of the total mixture is also 
obtained in the divergence form 



d_ n 

dt 



Qpu 2 + e - a w 2 ^ + div ^Qpu 2 - a w.u - aj u - U + e uj = 0. 
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, _ _ _ dp d\p A . . , , . 

with U = — $ H — $1 + — V and it is possible by a similar calculation 

ere at at 

than in section 3 to deduce an acceleration wave associated with the entropy 

gradient. 
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